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C^ ■ Abstract 

CIhI We discuss the finite size behaviour in the canonical ensemble of 

^ ■ the balls in boxes model. We compare theoretical predictions and 

numerical results for the finite size scaling of cumulants of the energy 
distribution in the canonical ensemble and perform a detailed analysis 
^ ' of the first and third order phase transitions which appear for different 

j^ . parameter values in the model. 

Introduction 

The balls in boxes model is a simple statistical model describing an ensem- 
ble of balls distributed in boxes subject only to a single global constraint. 
In some circumstances the model undergoes a phase transition driven by a 
condensation of the balls into a single box. 

The model has been introduced as a mean field approximation to lattice 
gravity |[|, 0] . Despite its simplicity, the model captures the main features of 
the phase transition observed in lattice Euclidean quantum gravity models 
[ie dynamical triangulations) such as the discontinuity of the transition |Q, ^ 
and the appearance of singular structures |^, H, |^, || . 
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Variations of the phase transition of the balls in boxes type can be found 
in many areas of physics like the thermodynamics of hadrons or strings |P, |10| 



branched polymers [^, |T^ or percolation]^. The balls in boxes model is also 



closely related to urn models [T^, |14[ and the spherical model |T^ . It is thus 
a rather generic type of phase transition. 

The model has an interesting phase structure : by tuning a single pa- 
rameter one can change the order of the transition or make the transition 
disappear entirely [^ |T^. This sort of phase transition was discovered in the 
spherical model by Berlin and Kac |l^ . The nature of the transition is quite 



different from more standard phase transitions in the theory of critical phe- 
nomena which result from the appearance of long range spatial correlations. 
Here the transition does not refer to any correlations in space. Instead, it is a 
kinematic condensation which comes about when a fugacity of the series rep- 
resenting the partition function hits the radius of convergence of the series. 
This fugacity sticks to the radius of convergence and cannot move. At this 
point the corresponding physical system enters the condensed phase. From 
this point of view the mechanism of the transition is very similar to Bose- 
Einstein condensation. The basic difference between the balls in boxes type 
of condensation and the Bose-Einstein one is that in the latter the system 
condenses into an energetically favoured ground state, while in the former 
balls condense into a box which is indistinguishable from the remaining ones. 
The balls in box condensation thus spontaneously breaks the permutation 
symmetry of the boxes. 

The thermodynamic limit of the balls in boxes model and the relation 
to lattice gravity has been discussed in a series of papers [|l], ||, |16[. Here 
we discuss finite size effects and confront theoretical predictions with the 
numerical analysis. The numerical analysis is done using an exact algorithm 
which recursively generates partition functions for systems of moderate size 
ie up to a few thousand boxes. 

The model and thermodynamic limit 

In this section we briefly recall the model P, |l^. The model describes N 



balls distributed in M boxes. Each box has a certain weight p{q) which 
depends only on the number g, of balls in it. The weights p{q) are the same 
for all boxes. For a given partition of balls in boxes {gi, g2, • • • , 5'Af }, the total 
weight of the configuration is equal to the product p(gi) . . -piqAi) - ie the 

"'^The appearance of a percolating cluster in the percolation theory can be translated 
into the condensation of balls in one box, which after condensation contains extensive 
number of balls. 



boxes are almost independent. They are not entirely independent because 
of the constraint qi + ■ ■ ■ + qm = N. It is this constraint which makes the 
model nontrivial. The model is defined as a sum over all partitions of A^ 
balls in M boxes weighted by the product weight. This is a fixed density 
ensemble, with density p = N/M. In the thermodynamic limit M — *■ cxo, 
and p = const, the behaviour of the system depends on the density p. For 
some choices of weights p{q), there exists a critical value pcr above which the 
condensation of balls takes place, which means that a box appears in the 
system which contains a finite fraction of the balls. The phase transition in 
p corresponding to this condensation is of third or higher order . 

The fixed p ensemble is a kind of microcanonical ensemble. One can also 
consider ensembles with varying density. There are two natural candidates 
for such "canonical" ensembles : for a fixed M one can allow for varying N by 
introducing a conjugate coupling to A^. In this formulation boxes decouple 
and we are left with M copies of the urn model. Alternatively for fixed A^ 
one can allow for varying M. This corresponds to the ensembles studied in 
lattice gravity [0] and will be the subject of our investigations. 

The partition function for this canonical ensemble reads : 

oo 

M=i qi,---,qM 

oo 

= ^ e''^^2(Ar, M) . (2) 

M=l 

The canonical function can formally be treated as a discrete Laplace trans- 
form of the microcanonical partition function z{N, M) as the second line of 
the equation above shows. Analogously, in the same language, we can intro- 
duce the grand-canonical partition function as the Laplace transform of the 
canonical one : 

oo 

Zip, ^) = Y, e-^'^ZiN, k) . (3) 

N=l 

Each box is assumed to contain at least one ball. Therefore the sum in the 
last equation starts from A^ = 1. 

The weights p{q) are a priori independent parameters that we assume 
to be non-negative. In order that the large A^ limit be well defined the 
weights cannot grow faster than exponentially with q. On the other hand the 
exponential growth factor is irrelevant from the point of view of the critical 
behaviour of the model, as can be seen below. Namely, if one changes weights 



as : 

p{q) -^ p'{q) = e-'^°e''«''p(g) , (4) 

where /iq and Kq are some constants, the partition functions change as : 

z{N, M) -^ z'{N, M) = e''°^'''''^'^z{N, M) 
Z{N, k) ^ Z'{N, k) = e''°^Z{N, k - kq) (5) 

2:(/i, k) -> Z'(/x, k) = Z(/x - /io, K - Ko) • 

Thus the change of the exponential factor in the weights merely causes a 
redefinition of the coupling constants /z and k. So from the point of view of 
the phase structure of the model, only the sub-exponential factors of weights 
matter. Here we restrict ourselves to the power like weights : 

p{q)=q-P, g = l,2,..., (6) 

which capture a variety of essential phase transitions in the model. Some 



remarks on the behaviour for weights of general form can be found in |T6 . 

In the large A^ limit the model can be solved analytically [|16|. The free 
energy : 

<P{k)= \im i-logZ(iV,«) (7) 

has a singularity in n at Hcr = —^og({P), for /? G (1, oo) and has no singu- 
larity otherwise. ({(3) is the Riemann Zeta function. The phase for k > n^r 
is called fluid, and for k < K^r condensed. 

The order parameter for the transition is the first derivative of the free 
energy with respect to k which corresponds to the inverse average density of 
balls per box : 

It vanishes in the condensed phase ie for k < /£„•• The average (...) is taken 
with respect to the canonical ensemble @. 

The order of the transition depends on (3. For (3 in the range (3 G (2, cxd) 
the transition is first order. When n -^ «;+ approaches the critical point 
the order parameter goes to a nonzero constant raise = C{P)/({P ~ 1) corre- 
sponding to the discontinuity of the order parameter at the critical point. If 
one treats M in the partition function (^ as the energy of the system, then 
the discontinuity of r corresponds to the latent heat. For /3 — > 2 the latent 



heat vanishes and the transition becomes continuous. For (3 in the range 
/3 G (1, 2], r has a branch point singularity at /€„•• In this case, the order pa- 
rameter vanishes when one approaches the critical point Ak = k — n^r -^ 0"^ : 

2 — 3 
r = 0(^) ~ Ak^^, where Xi = — . (9) 

The transition is ra-th order for /5 G ("^^j^^ri]- When /3 approaches one 
(/3 —* 1+) the transition becomes softer and softer in -^ oo) and eventually 
disappears at /3 = 1. For /5 smaller than one the system has only a fluid 
phase and displays no phase transition. 



Finite size scaling 

Let us denote the most singular term in the free energy by 4>sing{i^) and the 
related exponent by xq : 

<Psing ~ A/€^° . (10) 

To match it with the scaling (^ for /9 G (1, 2), we must have : 

.0 = ^. (11) 

For large but finite N one expects that this behaviour gets substituted 
by the double scaling law with two exponents which we denote by Aq and 



i?(see [jT8[ and reference therein) : 

(f)sing{AK, N) = iV^V(A«:iV^) , (12) 

where /(^) is a certain universal function of one argument ^ = AkN^ . When 
Ak is fixed and A^ goes to infinity one expects to asymptotically recover the 
singularity of the thermodynamic limit ^^^^(Ak, A^ -^ oo) ~ Ak^° indepen- 
dent of N. This means that the universal function must have the following 
large ^ behaviour : /(O ~ C^° and : 

Xo = -Ao/B. (13) 

The singular part of the free energy is mixed with less singular and an- 
alytic parts. While for An — > the (psing will eventually dominate over less 
singular parts it can still itself be dominated by the analytic part. It is there- 
fore difficult to extract the scaling part 4>sing directly from the function 0. 
Instead, one does it by calculating the n-th derivative of 4> for n large enough 



that the exponent Xn is negative. We use the convention that the exponent 
Xn is the power <p'2nq{'^) ~ Ak^" of the most singular part of the n-th deriva- 
tive. Thus we trivially have a;„ = Xq — n. The first value of ra = 1, 2 . . . for 
which Xn is negative, gives the order of the transition. For this n, the singular 
part of the n-th derivative blows up at Kcr and dominates the analytic part. 
Therefore in this case one can skip the subscript sing when writing a scaling 



formula analogous to (|T2|) for divergent derivatives : 



0(")(Afi;,Ar) = Ar^"/(")(A/tA^^) + ... , with A^ = nB + Aq (14) 

keeping in mind that there are corrections, (denoted by dots) which may 
be important for finite A^. Asymptotically, for N —>■ oo they are, however, 
negligible in comparison to the displayed part. 

The knowledge of the exponent Xq does not suffice to calculate Aq and B. 
In the standard finite size scaling analysis the value of the exponent B can 
be obtained by a simple argument |T^, |18[. There are two relevant length 
scales in a d-dimensional system undergoing a continuous phase transition : 
the linear extension L = N^/'^ and the correlation length ^, diverging as 
^ ~ (AT)^", when AT = T — T^r goes to zero. The critical behaviour sets 
in when the lengths become comparable : 

L^i (15) 

giving : 

AT ~ N'^ — > B = — . (16) 

Alternatively one can say that critical properties of the system depend on 
the dimensionless ratio : ^/ L. We use here a similar reasoning. Namely, 
we extract characteristic scales from the distribution of the box occupancy 
number 7r(g) ie the probability that a box has q balls : 

In the large A^ limit the distribution has the form : 

7r(g) ~ p{q) g-'^^^'^)" (18) 

as can be seen from (|^. There are two scales which govern this distribution 
when the critical point is approached from the fiuid phase : the damping 
scale factor [(J){/S.k)Y^ and the system size A^. Thus, in analogy to ( [T^ ) one 



expects that the dimensionless combination of the scales which defines uni- 
versal critical properties of the system is : [(j){An)]^^/N. Using the relation 
dmi ) we eventually get : 

B = Xo^ = P-1 and Aq = -1 . (19) 

Inserting this into (0) we obtain An = nB — 1. This is a sort of Fisher 
scaling relation^. 

There are some thermodynamic inequalities which follow directly from 
the Fisher relation. Firstly, one can see that B must be greater than zero, 
-B > 0, in order that there exist n such that A„ > 0. Otherwise there were 
no transition. This inequality is in accordance with the inequality (3 > 1 
which gives the condition for existence of the transition in this model, as 
discussed in the text after the equation (|^). Secondly, A2 = 2B — 1 must 
be smaller than one A2 < 1 which amounts to -B < 1. This is because the 
second cumulant may not grow faster than ^iV as can be seen from (^2|) by 
taking into account the fact that values of M lie in the range [1, iV]- In fact 
the limiting value B = 1 corresponds to the first order transition scaling 
which follows from the presence of a non-vanishing latent heatQ ie rdisc > 0. 
This holds for all /3 < 2. In this range of (3 the exponents A2 = B = 1 are 
constant, while the discontinuity r^isc changes with /?. 

Formula (0) gives a practical way of computing the most relevant sin- 
gularity of the partition function from the finite size analysis of numerical 
data. Namely, one computes divergent derivatives, for them estimates the 
exponents An, B and using the relation ^4^ = ^10 + nB one calculates : 

Xn = -An/B = xo + n. (20) 

Numerically, derivatives 0*^"-' of the free energy are computed as cumulants 
of the canonical distribution for M (B) : 



^(1) 



4- - <^^y<^ (22) 

^,3, _ (A/''>-3(AO(M) + 2(M)' ^ pg^ 



^In the standard considerations of the second order phase transition in the theory of 
critical phenomena the exponent B ~ l/ud and A2 — a/vd^ The exponent a is the heat 
capacity exponent ie the second cumulant of energy. Thus, the relation A2 — 2B — 1 
corresponds to a = 2 — lyd, known as Josephson's or Fisher's law. 

"^In the standard theory of critical phenomena B = l/ud = 1 gives the canonical 
exponent ly = 1/d. 
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Figure 1: The second and third derivatives of (f){K,,N) at /5 = 1.4 in ther- 
modynaniical hniit and for two different system sizes. The critical value is 
at Kcr ~ —1.133, the position of the maximum of the second cumulant is at 
i^max ~ —0.865, and the minimum of the third one at K^m ~ —0.449. Those 
are indicated by dashed vertical lines. 



For finite A^ the first cumulant (p^^^K, N) grows monotonically from zero to 
one as a function of k, while the higher order cumulants 0'-"^ have n — 1 ex- 
trema. If a cumulant is divergent, the leftmost maximum corresponds to the 
scaling, critical part described by the formula (|T3p. The remaining extrema 
lie in the fluid phase and come from the non-scaling part. In the thermody- 
namic limit the leftmost maximum of a divergent cumulant approaches the 
singularity : (J)^'^\k) ~ Ak^". 

As an example consider the figure |l| where derivatives of (f){K, N) are plot- 
ted for /3 = 1.4. In this case system undergoes third order phase transition at 
Kcr = — log C(1.4). When k, approaches the critical value An = k — Hcr -^ O"*" 
the first and second derivatives of 0(k) vanish as 0^^-* ~ A/t^/^, 0^^-' ~ Ak^^"^, 
whereas the third one diverges as 0^^-* ~ A/t^^/^. The second derivative has 
a maximum that lies in the fluid phase far from the transition (|T]). The third 
derivative has a minimum in the fluid phase away from the critical value Kcr- 
For flnite size the third derivative (p^^^n, N) has a maximum and a mini- 
mum. When N goes to inflnity, the position of the maximum tends to Kcr 
and its height grows according to a flnite size scaling ([T^). The position of 
the minimum approaches a value Kmin which lies in the fluid phase far from 
the critical region. 
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Figure 2: The probability distributions of the energy r = M/N for pseudo- 
critical values of K^r{N) at /? = 2.6 for A^ = 256, 512, 1024, 2048. 



For higher order cumulants the number of such noncritical extrema in- 
creases with the order of the derivative. It may even happen that for a given 
system size the height of a non critical maximum is larger than of the critical 
one. Generally, to determine A^ and B one should analyze only the scaling 
of the leftmost maximum of a cumulant, where the information about the 
singularity of 0^"^ is encoded. 

In the next sections we will present a numerical analysis of the finite size 
data for the different ranges of (3. Before doing this, let us briefly describe 
the algorithm to generate the finite size partition function is Z{N, k) 

We calculate the partition function Z{N, k) in two steps. First we com- 
pute values of the microcanonical partition function z{N, M) for all M in 
the range [1, A^] by the following recurrence relation: 



N 



z{N,M) = J2^iN -q,M -l)p{q), 



(24) 



with the initial condition ^(l,^) = q~^. Inserting the numbers z{N,M) 
directly to the definition @ we obtain Z{N, k). The maximal size N which 
can be reached by this procedure is a few thousand. It is limited by floating 
point instabilities accumulated in the recurrence relation (^41). To test the 
stability we check whether the results stay intact (up to a shift in n) under 
the change of weights (H). 




Figure 3: The probability distributions of the energy r = M/N for pseudo- 
critical values of K^riN) for A^ = 2048 at /3 = 2.2, 2, 4, 2, 6, 2.8, 3.0. 

Finite size analysis results 

First order phase transition 

We start the finite size analysis with the range (3 > 2. According to the 
discussion above, the phase transition should be first order. A typical signal 
of first order phase transition is a double peak in the distribution of the 
energy corresponding to the coexistence of two phases at the transition with 
the latent heat being the separation between peaks. Let us look for this 
signal in our case. In figure || we plot the pseudocritical distributions of the 
energy density r = M/N for different system sizes A^ for /3 = 2.6 for values 
of K,cr{N) for which the both peaks have equal heights. These values can be 
taken as pseudocritical ones. The position of the left peak goes to zero when 
A^ goes to infinity, while of the right one to rdisdP) = C{P)/({P ~ 1) = 0.571 
known from the analytic calculations in the thermodynamic limit. When /3 
changes the position of the right peak moves (figure |^). For large (3, the 
discontinuity raise goes to one, while for /3 — >• 2 it disappears and the two 
peaks merge. The depth of the valley between the peaks increases with A^ 
which means that the configurations which do not belong to either of phases 
become more and more suppressed. The suppression becomes more visible 
for larger (3 (figure ^. On the contrary, for /3 close to 2 the valley is small or 
even absent . In this case the size of the system must be increased sufficiently 
for the valley to be visible. 
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Figure 4: The pseudocritical value KcriN) versus 1/A^ for (3 = 3.0. The 
curve going through the data points corresponds to the best fit k = n^r + 

a/N (l + b/N). 

The range of k for which two peaks coexist is called pseudocritical region. 
For K in this region the relative heights of the peaks vary. The extent of 
the region is inversely proportional to the size of the system, A^. Outside 
the pseudocritical region only one peak survives. If k, changes within the 
pseudocritical region, the average of the distribution r = {M)/N moves very 
quickly between two peaks. The slope of the curve r{K,) grows linearly with 
A^, and eventually becomes infinite when A^ is sent to infinity, leading to 
the discontinuity raise at Hcr- The slope of the curve corresponds to the 
second derivative, which is heat capacity. Because the heat capacity grows 
linearly with N, the system has latent heat corresponding to the energy 
needed to move states of the system from one peak to the other. Another 
characteristic signal of first order phase transition which can be read off 



from the formula (14) is that the position of the maximum of the second 
cumulant should asymptotically lie on a curve AnN = const, which means 
that k{N) = Kcr + const/N. As an example we show in figure ^the behaviour 
of Kcr{N) versus A^. The data points are fitted to the formula with next-to- 
leading corrections of the standard form : k(N) = Kcr + j^i^ + j^)- The fit 
give^ Kcr = -0.185(1), a = -15.9(3) and b = -24(2). The coefficient b of 



"^To estimate the errors of the fit parameters we use the following procedure. The 
data consists of n points for different volumes. We successively omit one of them and fit 
the formula to the remaining ones obtaining n different fits. Having done this, for each 
parameter of the fit we have a distribution of n values. We find the average and the 



11 



the correction term is large. Skipping the correction h/N in the fit would 
decrease the quality of the fit and would significantly change the estimate of 
the Kcr which agrees with the analytic result K^r = — logC(3) = —0.184 with 
the correction term. 

To summarize, the standard signals and finite size scaling characteristic 
of a first order transition are observed in the range /3 G (2, cx)) as expected. 

Continuous phase transition 

For (3 G (3/2, 2) the linear growth of the maximum of the second cumulant 
of the first order phase transition changes to a sublinear behaviour ~ N"^'^ , 
where < yl2 < 1, corresponding to a second order phase transition. There 
is no double peak signal in the distribution of M. The pseudocritical point is 
defined by the value of n at which the second cumulant of the distribution of 
M is the largest. The second order phase transitions have been extensively 
analyzed numerically in a number of papers. Therefore, we prefer here to 
go directly to the third order phase transition, (3 G (4/3,3/2), where some 
new ingredients like the presents of non-scaling extrema appear on top of 
the standard finite size effects known from second order transitions. As 
an example, we consider the case /3 = 1.4 mentioned before, for which we 
expect that the second cumulant does not diverge, while the third one does 
with the exponent X3 = —1/2. As discussed before, we will concentrate 
our attention on the leftmost maximum of the third cumulant, which signals 
the appearance of the divergence in the thermodynamic limit. For A^ going 
to infinity, the height of the maximum is expected to grow as a power A^'^^ 
(|T^), possibly with some deviations for finite A^ coming from non-scaling part 
denoted by dots in (0). Indeed, it turns out that in the range of A^ from 16 
to 4096 which we have covered by the recursive method (^^ the corrections 
to the asymptotic formula are strong, as can be seen in figure |^, where the 
height of the maximum of the third cumulant versus A^ in the log-log scale 
is plotted. This is clearly non-linear and we define the effective exponent 
A'i{N) as the slope of the line fitted through three consecutive points with A^ 
standing for the biggest size A^ of the three points used in the linear fit. The 
effective exponent is clearly far from reaching an asymptotic A^-independent 
value (figure |). To be able to find the asymptotic value A*^ = A-^i^N ^>- cxd) 
from this data, one would have to know the form of subleading terms. Since 
we do not, we propose here a phenomenological approach. We postulate a 



width of the distributions which we take as the mean and the error of the parameters. In 
doing this we assume that aU the data points are equally important, which is the simplest 
possible assumption. This procedure is not a statistical analysis, but merely a way of 
presenting data. 
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Figure 5: The data points represent the maximum of the third cumulant 
for different A^. They are plotted in a log-log scale. 

phenomenological form of corrections to the running effective exponent 

A^iN) = A; + bN-". (25) 

It turns out, that this formula gives a good fit (figure P) with the values 
Al = 0.16(1), b = 1.58(1) and c = —0.27(1), estimated by the same procedure 
as discussed in the footnote on one of earlier pages. From the fit one can make 
a qualitative estimate of the subleading corrections. A simple calculation 
shows that A^ must be of order of 10^ in order that the running exponent 
differ from the asymptotic value Ag, say, by 0.04. It is of course a rough 
estimate, but it points to how slowly the corrections change with A^. 

The standard method of measuring the exponent B relies on tracing the 
position of the maximum for different sizes. For large A^ it approaches an 
asymptotic value. This value corresponds to the critical temperature. The 
exponent in the function which measures the distance of the position of 
the maximum from the critical temperature is proportional to An ~ N~^ . 
This gives the possibility of computng B. Due to the strong subleading 
corrections, this method does not give a good estimate in our case. Instead 
we propose another one. Unlike in the standard Monte Carlo, in our case the 
cost of computing cumulants does not grow with the order of the cumulant. 
We use this fact, to compute the growth exponents An for higher order - 
divergent cumulants, n = 3, 4 . . . . Then we use the formula for linear growth : 
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Figure 6: The effective exponent A^lN) calculated as a slope of the lines in 
the previous figure. 

An = Aq+uB, which allows us to compute the exponent B as the slope of the 
line. As an input for An we take the asymptotic values ^* with the errors, 
for n = 3, 4, 5, 6 estimated by the method discussed above. The fit is shown 
in figure |. It gives : Aq = -1.005(42), B = 0.387(12) and x'^/dof = 0.45. 
The errors quoted correspond to the 99.7% confidence levelQ The deviation 
from the theoretical values Aq = —1.0 and B = 0.4 predicted from the Fisher 
relation are rather small : less than 1% for Aq and 4% for B . The former 
value agrees within the errors with the theoretical one, while the latter one 
is a little outside the error bars. 

However, when one inserts those values to the formula x^ = —Aq/B — 
3 and assumes the maximal correlation between Aq and B, the errors are 
strongly enhanced and one obtains X3 = —0.403(40). This is 20% off from the 
theoretical value X3 = —0.5 and twice as much beyond the estimated errors. 
The assumption about the strong correlation is dictated by the covariance 
matrix for the fit. On the other hand, if we instead assume that the Fisher 
relation indeed holds, as suggested by both the analytical and experimental 
evidence, we can set Aq = —1 and in this case we have only one parameter 
relation X3 = 1/B — 3 which leads to X3 = —0.416(80). 

In interpreting the results one has to remember that they depend on 



^ For a cross-check we have also done a fit assuming equal errors on each An and using 
condition C^/do/ = 1. The results Aq = -1.017(48) and B = 0.389(10) are consistent 
with the ones obtained before. 
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Figure 7: The linear fit of the exponents A^, 
the fine corresponds to the exponent B. 



for n = 3,4, 5, 6. The slope of 



the phenomenological fit (^3p extrapolated over several orders of magnitude. 
Such a fit gives meaningful results only due to the fact, that the small volume 
data points are exact, but on the other hand one has to keep in mind that in 
general such a procedure may introduce a systematic error which is hard to 
quantify. 

Reduction of the error would require either a knowledge of the form of the 
next-to-leading correction or pushing the computations to the system sizes 
of a few orders of magnitude larger. 



Discussion 

The balls in boxes model provides a useful laboratory for testing finite size 
scaling. On one hand one can predict critical exponents theoretically, on the 
other hand one can exactly calculate the finite size partition function. In this 
respect the model is exceptional, since most of the models must rely on noisy 
Monte Carlo data. A comparison of the asymptotic form (jT^) and the finite 
size results indicates the presence of strong corrections to scaling. The correc- 
tions come from the analytic and less singular part of the partition function. 
The corrections seem to be a very slowly varying function of the system size 
as figure ^ demonstrates. The phenomenological approach employed to es- 
timate the next-to-leading corrections and the limiting asymptotic value of 
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the effective exponent is possible because unlike the Monte Carlo data, our 
data points have no statistical uncertainties. 

The finite size analysis results support the hypothesis that the critical 
properties of the model are encoded in the one box probability vr(g). The 
universal scaling properties are obtained from comparing two scales present 
in the model : the inverse of an exponential fall-off parameter of the effective 
probability distribution 71(5) and the system size A^. The former one mea- 
sures the extent of fluctuations of the number balls in a box and the latter 
one serves as a natural cut-off of those fluctuations. As mentioned in the 
introduction this type of transition plays an important role in many physical 
systems. In particular it is relevant for understanding the transition in lattice 
gravity models^ 

The model provides a nice example of the finite size pattern for a first 
order phase transition. The results for the distribution are exact and easily 
reproducible so one can measure all relevant quantities and quickly test any 
new ideas. Thus, we hope, the model may also prove useful for studying 
generic properties of first order transitions. 
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